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Abstract 

Quantum optical phenomena are explored in artificial àtoms well known as semiconductor quan- 
tum dots, in the presence of excitons and biexcitons. The analytical results are obtained using 
the conventional time-dependent perturbation technique. Numerical estimations are made for a 
realistic sample of CdS quantum dots in a high-Q cavity. Quantum optical phenomena such as 
quantum Rabi oscillations, photon statistics and col·lapse and revival of population inversion in 
exciton and biexciton states are observed. In the presence of biexcitons the col·lapse and revival 
phenomenon becomes faster due to the strong coupling of biexciton with cavity field. 
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I. INTRODUCTION 



Recent progress in crystal growth techniques for the fabrication of nanostructure ma- 
terials has enabled us to realize the possibility of high-quality epitaxial layers enough to 
demonstrate excitonic optical properties of quantum dots (QDs) in the ultraviolet [1] to 
infrared [2] range including room-temperature excitonic lasing, high-temperature excitonic 
stimulated emission and high characteristic temperature for optical threshold power [3, 4]. 
Such nonlinear optical effects would be much more enhanced if biexcitons were involved 
in the optical processes because of the giant oscillator strength effect [5, 6]. In fact, low- 
threshold lasing based on optical processes associated with biexcitons has been theoretically 
predicted [7, 8] and observed [4]. 

Since its proposal by Purcell [9] and early work involving àtoms in cavities [10], cavity 
quantum electrodynamics has been actively pursued for its potential insight into fundamen- 
tal problems in light-matter interaction. Here, the energies and modes of photon states, to 
which electrònic transitions of àtoms can couple, are changed by the cavity confinement. 
Recently, efforts have been made to extend this work to the optical regime, which requires 
cavities with sizes of the order of few micrometers. This work has been made possible by 
the epitaxial fabrication of high quality layered semiconductors from which optical cavities 
have been made that have strong photonic confinement in one direction [11-13]. However, 
the preparation of optical microcavities which confine the propagation of light in all three 
dimensions is presently a challenge for modern microtechnologies. New techniques are be- 
ing developed to overcome the difliculties in developing such microcavities using quantum 
dot structures. Quasi atòmic light emitters in these resonators can be mimicked by semi- 
conductor quantum dots that exhibit discrete density of electrònic states. Arakawa et al, 
[4] fabricated QD làsers with vertically stacked self-assembled dots with a columnar shape 
which achieved a room-temperature threshold current of 5.6 mA, comparable to or even 
less than that of quantum-well làsers. The threshold current density has also been lowered 
to 90 A/cm 2 [12]. Pulsed and continuous-wave operation of 1.3 fim quantum dot làsers at 
room temperature has been realized, which is a milestone towards the application of QDs 
to ultrafast fiber òptic data transmission and optical interconnects [2]. Artemyev et al, [14] 
have developed high-Q microcavities using tiny hollow microspheres made of the transpar- 
ent polymer polymethylmethacrylate (PMMA) with smooth, highly reflecting surfaces that 
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exhibit all properties of an optimum optical cavity. Such improved làser performance has 
now enabled physicists to study various unique properties of interaction between QDs and 
light. 

The study of interaction of quantized radiation with matter is important in QDs, be- 
cause the typical sizes of QDs are much smaller than the wavelength of light being used 
(for example the typical sizes of QDs are CdS - 1.7nm [15], and CdSe ~ 4nm [16]). The 
smaller cross sectional areas of these QDs allow only few photons to interact with each of 
them. Hence, the application of quantum optical tool to QD may yield interesting results. 
Accordingly, with the recent developments in the fabrication of high-Q microcavities hav- 
ing 3-dimensional photonic confinement on one hand and the enhanced excitonic nonlinear 
optical characteristics of the QDs in 3-dimensional electrònic confinement on the other, we 
have made an attempt to explore the quantum optical phenomena in these structures in the 
presence of both excitons and biexcitons. 



II. THEORETICAL CALCULATIONS 



Time-dependent perturbation technique under the interaction picture has been employed 
to study the interaction of a single mode of cavity radiation with a small QD embedded in a 
high-Q microcavity. The atomic-like energy-level structure existing in these QD nanostruc- 
tures allows the photo-induced electrònic transitions to take place between the ground (|0)) 
and exciton (]e)) states. For large excitation intensities, this simplistic two-level picture is 
modified by the creation of biexcitons. Accordingly, in the present paper, while studying 
the nonclassical phenomena in QDs, we have incorporated the exciton and biexciton (|6)) 
states via a three-level system (3-LS) as schematically represented in Fig. 1. The allowed 
electrònic transitions (sòlid arrows) take place between the sates |0) |e) and |e) \b). 
The parity violating two-photon transitions (dashed arrow) between |0) =f± \b) are ignored 
in the present calculations. 

For small quantum dots of size (R) smaller than the exciton Bohr radius (cíb), the energy 
gap between the vacuum and exciton states is defined as [17] 

^ oe = ^ + ^-(^) 2 , (1) 
where fkv g is the band gap energy and n m i is mth root of the Zth order Bessel function with 
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m and / corresponding to the ls, lp, ld, ... 2s, 2p, 2d, ... levels of the electrons and holes. 
m r is the reduced mass of an electron-hole pair. Due to the finite Coulombic attraction 
between two excitons, the energy of the biexciton is reduced by an amount equal to the 
binding energy (Ai?) and defined by [18] 

frw eb = hu oe - AE. (2) 

We have assumed that at time t — 0, all the electrons are in the ground state. For time 
í > 0, the interaction of cavity photons with the QD creates one and two electron-hole pairs 
(i.e., excitons and biexcitons). Accordingly, the total wavefunction of QD-radiation can be 
defined under dressed state representation as 



I* (í)> = c ,n+i(í)|0, n + 1) + c e,n(í)|e, n) + ^ c 6>n _i(í)|6, n - 1). (3) 

ml ml 

Co,n+i(t), c etn (t) and Cb jn _i(í) are the probability amplitudes of the vacuum, exciton and 
biexciton states, respectively While writing the total wavefunction, we assume for the field 
mode that the initial state is coherent, such that 

\Cn(0)\ 2 = ^Ç, (4) 

where, n is the initial photon number. 

The interaction Hamiltonian of the exciton and biexciton states is represented by 

TLi = Tie + Tib-, (5a) 

with 

H e = -g oe [an eo e- iA ^ + a^ oe e lA ^] (5b) 

and 

H b = -g eb [an be e- lA ^ + ah eb é A ^}, (5c) 

where and a are the creation and annihilation operators of the electromagnètic field, 7Tjj 
= The detuning parameter is defined as A oe ( eb ) = uj — u) oe ( eb y g oe = \/i e E /h\y/n + 1 

and g eb = \/i b E/h\y/n are the coupling parameters of the exciton and biexciton, respectively 
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with the cavity field. Also E is the field per photon and 

li e = (e\fi\0) = n oe J <j)(r e ,r)dr (6a) 
fx b = {b\fi\e) = -VZfioe (6b) 

X 1 1 1 0( re > r > rh > r ï*^ re > rh ï drdr ' edrh - ( 6c ) 

for the same spin states. [19] Here, /i oe = [e\p oe \/(m ui oe )] is the transition dipole moment 
of the vacuum to exciton state transition and m is the free electron mass. 0(r e ,r) and 
<p(r e ,r,r h ,r) are the exciton and biexciton wave functions, respectively. They are related 
to the single particle wave function 0(r) as [18, 19] <f>(r e ,r) = <f>(r e )(j>(r) and <p(r e ,r,r h ,r) = 
<f)(r e )(f)(r)(f)(r h )(f)(r), with 

ó(r)= F^S rtKmr/R) 
nn V R* j l+1 (K nl ) ■ 

Using Schrodinger's equation under the interaction picture, the equations of motion of 
probability amplitudes are found to be 

ào,n+i(t) =i^5f oe c ein (t)exp(iA oe í), (7a) 

e 

Ce,„(í) =Í9oeCQ,n+l(t) exp(-íA oe í) 

+ igebCb,n-i(t)exp(iA eb t), (7b) 

C b ,n-l(t) =ÍgebC e , n (t) exp(-iA eb t). (7c) 

The above set of coupled equations can be solved by assuming 

c o ,n+i(0 = Aex.p(iÇlt), (8a) 
Ce, n (í) = B exp[i(ü - A oe )t], (8b) 
C6,n-i(í) = Cexp[i(n - A ob )t], (8c) 

where A, B and C are the time-independent constants, Vt is an unknown parameter having 
the dimension of frequency and A ob = A oe + A eb . Use of (7) and (8) yields a cúbic equation 
in fl as 

íl 3 + a 2 ü 2 + ajn + al = (9) 
with al = -|5f oe | 2 A o6 ; a\ = A oe A ob - \g oe \ 2 - \g eb \ 2 and a 2 = ~(A oe + A ob ). The solution of 



5 



íl are found to be 



fi 2 

3 



a 2 X 2 - 12Y 2 



02 
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X 2 e^ 



12Y 



6X 



(10a) 
(10b) 



with 



X 3 = 360102-10800-803+^^40? - a\a\ - \%a\a 2 a\ + Yïa% + ^a\a\ and Y 2 = a\--l. 

v O 

A. Rabi Oscillations 

For analytical coníirmation, we have reduced the generalized equations obtained for a 
three-level system (3-LS) to a conventional two-level atòmic system (2-LS). For a 2-LS, in 
the absence of biexcitons the coefficients are corrected by substituting g e b = and A & = 
and it is found that a = 0, a\ = — \g oe \ 2 and a 2 = — A oe , which reduces the cúbic 
equation to a quadratic equation in íl We find that íl 2 vanishes in the absence of biexcitons 
and the other two solutions remain finite. In a 2-LS, the magnitudes of these eigenvalues 
decide the freqüències of oscillation of the probability of the upper and lower levels (in 
the present case between vacuum and exciton states) [20]. This frequency of oscillation is 
widely termed as Rabi frequency. The three solutions of Í2 obtained in the present case of 
a 3-LS also exhibit a similar feature: ÍX,- obtainable from (10) may correspond to the Rabi 
oscillations occurring between the vacuum-exciton and exciton-biexciton states. However, 
in the presence of biexcitons, analytical results of Q could not be obtained due to the 
complex nature of the three solutions obtained in (10). Hence, we applied our analysis to a 
realistic sample of quantum dot of CdS embedded in an undamped microcavity. The material 
parameters of CdS taken from the experimental paper of Butty et al [15] are: huj g =2.56eV, 
AE = 28meV, R = 1.7nm, a B = 2.9nm and m h /m e = 4.2. Also, to make the present 
solutions compatible with the conventional eigenvalues [20], we define VL 1 - = ÍX,- — a 2 /3. The 
nature of dependence of Ví 1 - on exciton detuning parameter A oe are exhibited in Figs. 2 and 
3. In obtaining Fig. 2, we have incorporated the contribution of excitons only, while in 
Fig. 3, the contribution of both exciton and biexcitons are incorporated. The sòlid lines 
are obtained in the presence of finite cavity field and the dashed lines represent the same 
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obtained for n — 0. To increase the clarity of the curves, a thicker line has been drawn for fl 2 
while Çï[ and share same line-thickness. Parts of the dashed lines of Q[ (below resonance) 
and fl' 3 (above resonance) are not visible since fl' 2 for n = as well as n ^ overlaps with 
them. The curves obtained in the absence of biexcitons agree with the conventional results 
obtained for a 2-LS system.[20] For n — 0, all the curves vanish at A oe = 0. However, 
for n 7^ 0, the eigenvalues show repulsion at A oe = by an amount equal to ±g oe - The 
repulsive behavior of the two eigenvalues is attributed to the finite value of QD-cavity field 
coupling. The phenomenon of repulsion of the eigenvalues are also called as AC Stark shift 
or dynamic stark shift since the energy levels show shift due to the dynamic nature of the 
amplitude of the radiation. As expected from the analytical results for a 2-LS, one of the 
solution (Í2 2 = VL' 2 — a2/3) remains unchanged at zero. 

For n — 0, all the three curves obtained in the presence of both exciton and biexciton 
shown in Fig. 3, display biexcitonic signatures between A oe = and A eb = 0. A close 
observation of the curves obtained for n = yields the following: ^1,2,3 vanish neither at 
A oe = nor at A eb = 0. However, (i) Oi shows its minimum value at A oe + AE/2 and is 
repelled from the zero line by an amount AE/6; (ii) for A oe > < A eb , the second solution 
remains constant at —AE/3; for A oe < 0, it steadily increases to a value of AE/6 and again 
starts decreasing until A eb becomes equal to zero, and (iii) for decreasing detuning parameter, 
the final solution steadily increases to —AE/3 till A oe = 0; further, for A oe < > A eb it 
remains constant at —AE/3. For A eb < 0, it starts decreasing steadily. 

For n^0, dramàtic changes occur. Both as well as get repelled from Vt' 2 . In the 
presence of biexciton, Q! 2 plays a dominant role by repelling both Q[ and íl' 3 for n = and 
n^0. This leads to the conclusion that Q! 2 couples the biexciton states to the exciton states 
and hence the ground state. In addition to the shift along Y-axis, a small shift along X-axis 
can also be noted for all three curves in the presence of biexcitons. For n ^ the number 
of biexciton becomes finite, such that a\ as defined earlier modifies the resonance frequency. 
This leads to small red and blue shifts of the peak vàlues of the eigenvalues. These shifts 
indicate strong evidence for biexciton signatures in the eigenvalues Qi and Í2 2 . And 
does not show significant change in the presence and absence of biexcitons at A oe = and 
A eb = 0. 

The three solutions of Q obtained in (10) suggest the following modifications in the 
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assumed solution (8) 



Có,n+i(t) = ^ie íí!lí + A 2 e tU2t + A 3 e 



iQ3t\ —iA oe t 



c e ,„(í) = (5ie lS2lí + S 2 e lS22Í + B 3 e úht ) 
c b , n -i(t) = + C 2 e tÜ2t + C 3 e tn3t )e- iA ° bt . 



(lla) 
(11b) 
(11c) 



In order to calculate the nine time-independent unknown coeficients Ai, B { and Cj (i = 
1,2,3), we assume that, at t — 0, the QD is in the ground state |0) and the cavity field is 
in a coherent state. The corresponding boundary conditions at t — 0, are c 0)n +i(í) = c n (0) 
and c e , n (í) = C6 )fl _i(í) = 0. Use of these conditions in (7), (8) and (11) yields the solution of 
the unknown coeficients as 



( A 1 B 1 



V 



with M 



A 2 B 2 C 2 
A 3 B 3 C 3 

/i i i 

Çl\ ^3 



Cn(0) 

' \M\ 



^ 23 (a 2 3 + Ifi'oeH 9oeP23(^oe ~ I23) fioefie&A^ 

/3 3 l(«31 + Ifi'oel 2 ) fioeAsi(A oe - 731) g oe g eh (3 31 
^l 2 (ttl2 + Ifi'oeH fioe/3l 2 (A oe - 7i 2 ) fi oe fie6/9l2 ^ 



(12) 



«i 



íljíl,-, = ílj — ílj and 7^ = ílj + ÍX,-. For a 2-LS the 



^1 ^2 ^3 1 
above set of equations reduces to 



( Ai B 1 



A 2 B 2 C 2 
A 3 B 3 C 3 



( 



c„(0) 



\ 



2 fioe « 
^1 "fioe 




\ 



(13) 



In obtaining the above equations, we have assumed the biexciton-cavity coupling term (g eb ) 
to be zero. Comparison of (13) with the Standard solutions [20] for a 2-LS shows good 
agreement and confirms the validity of the model to a two-level system. 

Using (11)-(13), we have obtained the temporal variations of probabilities and exhibited 
these in Fig. 4. The sòlid, dashed and dotted lines represent |c 0jn _i| 2 , \c e>n \ 2 and |Q> irí+ i| 2 , 
respectively. The inset shown in the figure demonstrates the same parameters in the absence 
of biexciton contribution. In the absence of biexciton, the ground and exciton state popula- 
tions execute oscillation at a frequency equal to 2g e . Hence, the ground state is bleached for 
g e t = (2N + l)7r/2 with TV = 0, 1, 2, ... . As expected from earlier analytical simplifications, 
|o>,n+i| 2 remains constant at zero. In the presence of biexciton contribution, the boundary 



conditions are exactly followed. For t > 0, the ground state population starts decreasing un- 
til g e t = n/2. Simultaneously, the exciton population reaches its maximum when g e t = n/A, 
and later decreases to a minimum value at g e t = n/2. However, the biexciton population 
increases very slowly but attains its maximum value when g e t = n/2. This region is marked 
as 1 in the figure. This is an interesting region, since the population between |6) and |0) as 
well as between \b) and |e) is inverted. This population inverted region occurs again when 
g e t = 5n/2. The significance of the region 1 is that for g e t = (5N + l)n/2 and (N + l)5n/2, 
low threshold for exciton and biexciton lasing can be achieved. Another region of similar 
interest is marked as 2, where population inversion occurs between \b) and |e). All the three 
curves rephase for g b t = 4Nn. 

B. Photon Statistics 

The expressions (11) represent the probability amplitudes of the electron to occupy the 
states |0), |e) and \b) and the probability of n photons to occupy the states \n) at any 
arbitrary time t. The trace over the electrònic states yield the probability p(n,t) of the 
cavity field as 

p(n,t) = [Mt)\ 2 + |c e ,„(t)| 2 + |c M (í)| 2 ] . (14) 

ml 

Comparison of the above results with the semiclassical [21] results show disagreement, since 
the semiclassical techniques assume the photon number distribution at t = as c 0)ít+ i = 1. 
However, for an initially coherent state, the Poisson distribution allows the field to propagate 
in a coherent manner. In Fig. 5, variation of p(n, t) with time as well as number of photons 
is demonstrated for n = 10. Fig. 5a and 5b display contour plots in presence of excitons only 
and in presence of exciton as well as biexciton, respectively. The peak value is found to be 
at n. Both the curves strictly follow Poisson statistics with respect to n, but the envelope 
oscillates with time. A few peaks are obtained for smaller vàlues of n while the number 
of peaks starts increasing with n. The reason for these oscillations are attributed to the 
definition of g 0e which is proportional to \fn + 1. In the presence of excitonic contribution 
only, the photon statistics oscillates at a frequency g 0e . In the presence of biexcitons a clear 
case of chaos is observed as evident from Figs. 5c and 5d, which are obtained for n = 10. 
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C. Collapse and re vi val phenomena 



The phenomenon of collapse and revival of Rabi oscillations is studied in atòmic systems 
using atòmic inversion operator. The reason for the collapse and revival phenomenon is 
well understood now. The contributions corresponding to different n's interfere in such a 
manner that they initially go out of phase. After that they acquire a common phase, and this 
process is continuously repeated to obtain a series of col·lapses and revivals. In the present 
generalized model, we have assumed a three-level structure of excitons and biexcitons. The 
corresponding population inversion in these states can be found from 

W(t) 0e = ^(Kn+l(í)| 2 " |Ce, ít (í)| 2 ) (15a) 
e,n 

and 

W(t) eb = ^(|c e , n (í)| 2 - Iq^í)! 2 ). (15b) 

e,n 

In Fig. 6, the temporal variation of W oe is demonstrated in the presence (curve a) and 
absence (curve b) of biexciton contribution. Both the curves exhibit collapse and revival of 
population inversion between the ground and exciton states. A very slow revival is noted in 
the absence of biexcitons. In the presence of biexcitons, the strong coupling of the biexciton 
with cavity field reduces the revival time; hence a very small collapse period is observed. 
However, in the absence of biexcitons the QD behaves like a 2-LS; hence an atòmic like 
collapse and revival phenomenon is noted. In the absence of the biexciton contribution, W$ e 
oscillates around zero. With the biexcitonic contribution, it shows a positive DC shift. The 
transient behavior of excitonic and biexcitonic inversions, as exhibited in Fig. 7, does not 
show any remar kable change. In the presence of the biexciton Wo e , as well as W eb display 
the same collapse and revival time. However, W 0e oscillates around 0.3, while W eb oscillates 
around -0.2. 

III. CONCLUSIONS 

A simple quantum optical model has been developed to study nonclassical phenomena 
in semiconductor quantum dots in an undamped cavity. On reduction of the present cal- 
culations to a two-level system, these agree well with the conventional results. The results 
are applied to a realistic semiconductor QD of CdS of size 1.7nm. The transient nature of 
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photon statistics exhibits oscillations. The collapse and revival phenomenon in the presence 
of biexcitons becomes faster due to the strong biexciton-photon coupling. 
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FIGURE CAPTIONS 



FIG. 1 Schematic diagram of three-level structure of ground, exciton and biexciton states in a 
quantum dot. Sòlid lines represent allowed electrònic transitions, while the two-photon 
transition shown by the dashed line is not allowed. 

FIG. 2 The nature of dependence of eigenvalues Vt'j on exciton detuning parameter A 0e in CdS 
quantum dot in the absence of biexciton contribution. 

FIG. 3 Variation of Q'j with detuning parameter in the presence of biexciton contribution in 
a small quantum dot of CdS. 

FIG. 4 Temporal nature of probabilities of ground (dotted), exciton (dashed) and biexciton 
(sòlid) states in the presence of exciton and biexciton. The inset shows the same in 
the presence of exciton only. In the regions marked as 1, population inversion occurs 
between \b) and |e) as well as between \b) and |e). In the regions 2, inversion occurs 
between \b) and |e) states. 

FIG. 5 Photon statistics p(n, t) in a small quantum dot of CdS embedded in a microcavity for 
n = 10. The contour plot a is obtained in the absence of biexcitons while b is obtained 
in the presence of biexcitons. Plots (c) and (d) are obtained for n = 10. 

FIG. 6 The temporal variation of excitonic inversion (W oe ) in a QD of CdS. Curve a is obtained 
in the presence of biexcitons while curve b is obtained in the absence of biexcitons. 

FIG. 7 The nature of dependence of exciton and biexciton population inversion (W 0e and W e b) 
on time in a small quantum dot of CdS. 



13 



FIG. 1: Schematic diagram of three-level structure of ground, exciton and biexciton states in a 
quantum dot. Sòlid lines represent allowed electrònic transitions, while the two-photon transition 
shown by the dashed line is not allowed. 
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FIG. 2: The nature of dependence of eigenvalues Qj on exciton detuning parameter Ao e in CdS 
quantum dot in the absence of biexciton contribution. 
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FIG. 3: Variation of ílj with detuning parameter in the presence of biexciton contribution in a 
small quantum dot of CdS. 
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FIG. 4: Temporal nature of probabilities of ground (dotted), exciton (dashed) and biexciton (sòlid) 
states in the presence of exciton and biexciton. The inset shows the same in the presence of exciton 
only. In the regions marked as 1, population inversion occurs between \b) and |e) as well as between 
|6) and |e). In the regions 2, inversion occurs between \b) and \é) states. 
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FIG. 5: Photon statistics p(n, t) in a small quantum dot of CdS embedded in a microcavity for 
n = 10. The contour plot a is obtained in the absence of biexcitons while b is obtained in the 
presence of biexcitons. Plots (c) and (d) are obtained for n = 10. 
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FIG. 6: The temporal variation of excitonic inversion (W oe ) in a QD of CdS. Curve a is obtained 
in the presence of biexcitons while curve b is obtained in the absence of biexcitons. 
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FIG. 7: The nature of dependence of exciton and biexciton population inversion (Wo e and W^) 
on time in a small quantum dot of CdS. 
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